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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 2455 


A METHOD OF DESIGNING TURBOMACHINE BLADES WITH A DESIRABLE 
THICKNESS DISTRIBUTION FOR COMPRESSIBLE FLOW ALONG AN 
ARBITRARY STREAM FILAMENT OF REVOLUTION 
By Chung-Hua Wu and Curtis A. Brawn 


SUMMARY 

A rapid method for designing turbomachine blades of a given turning 
and a desirable blade-thickness distribution for a compressible non- 
viscous fluid flow along an arbitrary stream filament of revolution is 
presented. The method utilizes the guiding effects of the blade shape 
on the mean streamline shape and of the blade thickness on the specific 
mass flow along the mean streamline. After the flow on the mean stream¬ 
line is determined, the ‘extension of the solution from the mean stream¬ 
line to the blade surfaces is accomplished by the use of a power series. 
A number of blade profiles are obtainable for the total mass-flow 
requirement, and one is chosen for the best velocity distribution on 
the blade. The results obtained in the solution can be used for a 
direct check on the accuracy of series approximation and, also, for 
the more accurate determination of the velocity distribution along the 
leading and trailing edges of the blade. 

The method is illustrated with the design of several turbine cas¬ 
cades of highly cambered thick blades. The determination of the shape 
of the blades and the compressible flow past the blades was carried out 
by hand computation in 16 hours. One solution obtained by using three 
terms in the power series compared very well with an available direct 
solution and the blade circulation checks closely the specified turning 
angle. 

Because the surface of revolution, on which the blades are located,' 
is completely arbitrary, the method can be applied to axial-flow, 
radial-flow, and mixed-flow turbomachines. The variation in the normal 
distance between the stream surfaces of revolution can be taken into 
account, thus incorporating into the design the principal effect of 
three-dimensional flow. The method is readily applied to the design of 
channels on a plane and on a general surface of revolution. 
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IMRODUCTION 

The increasing use of compressors and turbines in aircraft power 
plants during the past 10 years has led to considerable research in the 
direct and inverse problems of two-dimensional potential flow past an 
infinite cascade of airfoils. In the inverse problem, the design of 
blades is often directed at control of the pressure or velocity distri¬ 
bution on the blade. Most of the methods are derived for axial-type 
turbomachines, in which the flow is assumed to take place on cylindrical 
surfaces (methods for incompressible and compressible flow are discussed 
in references 1 to 7 and 8 to 11, respectively). Methods for designing 
blades in a radial plane are given in references 12 and 10. 

In current axial- and radial-flow turbomachines, the flow surfaces 
are usually of a more general shape than either cylindrical or radial. 
Furthermore, the normal distance between adjacent flow surfaces varies 
along the flow path. A method was therefore developed at the Lewis 
laboratory for the design of blades for compressible flow along an 
arbitrary stream filament of revolution. Instead of the velocity 
distribution on the blade being the required result, the blade design 
is aimed at a desirable blade-thickness distribution required with 
respect to blade strength and Mach number in general and also the cool¬ 
ant passage requirement in the case of cooled turbine blades. The com¬ 
putation involved in this method is relatively simple and short, and 
the usual assumption of a linear pressure-volume relation for compres¬ 
sible flow is not required. 

For clarity and simplicity the method will first be given for 
compressible flow on a plane or a cylindrical surface and will be 
illustrated by a few examples. The method will then be given for the 
general case of compressible flow along an arbitrary stream filament of 
revolution. 


SYMBOLS 

The following symbols are used in this report: 

B differentiation coefficient 

H total enthalpy, h + ^ 

u 

1 2 

Hy. relative total enthalpy, h + -^ W 

Ct 0 

h static enthalpy 

v 
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_ , 1 tt 2 12 2 

I h + ^W -g 001 ’ 

L blade length projected on turbomachine axis 

l,op orthogonal coordinates on mean surface of revolution 

M mass flow 

P pitch or spacing 

p static pressure of gas 

r radial distance from axis of turbomachine 
S streamline 

t blade thickness in circumferential direction 
V absolute velocity of gas 

¥ velocity of gas relative to blade 

y distance in direction of pitch for plane flow and equal to 

for flow on cylindrical surface 

z distance along axis of turbomachine 

i W y - Wcp 

0 flow angle on stream surface. tan“ x -f- or tan“ x — 

w z Mi 

p density of gas 

a angle between tangent to meridional curve and axis, tan" 1 

T normal thickness of stream filament of revolution 

X ratio of specific heats 

\1 ! stream function 

0 ) angular velocity of blade 

Superscript: 

* 


dimensionless value 
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Subscripts; 
e exit 

i inlet 

Z,cp meridional and circumferential components 
m mean streamline 

p pressure surface of blade 

s suction surface of blade 

T total or stagnation state 

y y-component 

z z-component 


DESIGNING BLADES FOR COMPRESSIBLE FLOW IN PLANE 
OR ON CYLINDRICAL SURFACE 
General Description of Method 

In a recent investigation of compressible flow through a typical 
cascade of turbine blades (reference 13), the following results were 
obtained: 

{l) The shape of the mean streamline follows approximately that of 
the mean channel line of the cascade and has a lower curvature. 

(2) The variation of the ratio of the specific mass flow on the 
mean streamline to its inlet value follows the trend in the variation 
of the ratio of pitch to channel width (inside the channel the ratio of 
the specific mass flow is about 4 percent greater than the ratio of 
pitch to channel width). 

(3) The variation in fluid properties across the channel can be 
represented by a second-degree function for engineering accuracy. 

These results were used herein to develop a rapid method for 
designing cascade blades for either compressible or incompressible flow. 

This method of blade design starts with the calculation of the com¬ 
pressible flow on the mean streamline. With the flow angle upstream 
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and downstream of the blade determined by the velocity diagram and the 
pitch and the axial length of the blade given, a particular mean stream¬ 
line shape may be sketched by the designer. The designer may use the 
shape of the mean streamline itself as a parameter of the cascade or 
may specify the mean streamline shape according to figure 17 or refer¬ 
ence 13 (or even better according to any available relation of a simi¬ 
lar blading) to lead to a certain blade camber line. (if only a number 
of points of this streamline are specified, it is important that the 
values given are numerically smooth) . In addition to the mean streamline 
shape the designer further specifies at a finite number of points along 
the mean streamline (such as z^ . . . z^g in fig. l), the ratio of 

specific mass flow on the mean streamline to inlet value. The values 
of the specific mass flow on the mean streamline are determined by the 
blade-thickness distribution which is desirable from the consideration 
of blade stress and Mach number in general, and the consideration of 
the additional requirement of coolant passage in the case of cooled 
turbine blades, and by a relation between blade-thickness distribution 
and specific mass flow on the mean streamline, such as shown in fig¬ 
ure 19 of reference 13. 

With these specified values, the velocity components and the density 
are very easily determined at the specified points on the mean stream¬ 
line. The variation of velocity components and density in the pitch 
direction are then obtained by using power series in that direction. 

The derivatives in the series are determined from the fluid state on 
the mean streamline by the use of equations of continuity and motion 
and the density-velocity relation for isentropic flow. A number of 
blade profiles and their velocity distributions are obtained by inter¬ 
preting the starting mean streamline as dividing the inlet mass flow 
into two slightly different amounts in the channel. The velocity dis¬ 
tributions on the blades are compared and the best one is chosen. If 
the blade shape, the thickness distribution, or the velocity distri¬ 
bution around the blade obtained requires some modification, the 
values specified for the mean streamline should be modified and the 
process repeated. Because of the relatively short computation involved 
(in the illustrative examples of turbine blades, only 16 hours were 
required for the compressible solution by using the first three terms 
in the series and 6 to 10 stations inside the channel), modifications 
of the solution for more desirable blade thickness or velocity distri¬ 
bution is practical. Families of blade elements can be built up very 
quickly this way for any particular application. 


Basic Relations 

The steady two-dimensional isentropic flow of a nonviscous fluid 
in a plane or on a cylindrical surface is governed by the following 
equations of continuity and irrotational motion, and the isentropic 
pressure-density relation: 
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£)(PW Z ) d(pWy) 

5z f 5y ~ ® 

dw z dw v 

5T'S^'° 

P = Kp T 


( 1 ) 

( 2 ) 

(3) 


In these equations, the z coordinate is chosen along the machine axis 
and the y coordinate is chosen along the pitch direction (y = rep in 
the case of cylindrical flow with r equal to a constant)• 


Consider first the gas flow along a streamline somewhere in the 
midpart of the channel formed by two neighboring blades, such as ab 
in figure 1. The coordinates of the streamline and their differentials 
are related, respectively, by the following two equations: 


S(z,y) = 0 


(4) 


dS 


dz+||dy 


= 0 


(5) 


When the variation of the fluid state along the streamline is considered, 
any quantity q on the streamline is a function of z only, that is, 

q = q. £,y(z)] (6) 

The total derivative of q with respect to z is 


dq _ Sq bq dy 
dz 3y dz 


But 


as 

dy 3z 
dz ~ “ Ss “ 

5y 

Hence equation (7) may be written as 



tan p 


dq _ dq 
dz 3z 


tan p 


(7) 


( 8 ) 


O) 


When equation (9) is used, the continuity and irrotationality relations 
can be written as 
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a(pwj 


- tan fi 


B(pW z ) d(pW v ) 


dW v dW„ dw v 

IT - 5T ' tan ^ ' 0 < 

The variation of density throughout the flow region can he most 
conveniently expressed in terms of its inlet value through the use of 
equation (3) as follows: 


p* = = 

^i 


P± P T,i 


■2V-1 


A tabulated general relation of p/p^ ^ in equal intervals of W 2 /}^ 
can first he calculated, from which either a table for p/p. in equal 
intervals of W or a graph can be.easily constructed for each indi- 
vidual case and used for the evaluation of density from the velocity. 

Along the chosen (mean) streamline, where the slope is known and 
pW z is given at a number of stations, the density at these stations can 

be obtained as follows: Rewrite equation (12) as 


p* = £- = 
Pi 


1 - ^2 (P*V se c P) 

V 


(12a) 


where 


(p*W * sec p) 


1 - (P*) 


* (r-i) 


(12b) 
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Either equation (12a) is used to prepare a table of p * fcJr equal inter¬ 
vals of (p*W z * sec p) by an iterative process, or equation (l2b) is 
used to compute p*W * sec p for a number of values p*, which are 
then plotted as a graph. After densities have been obtained from either 
a table or graph, the velocity components on the chosen (mean) stream¬ 
line are readily computed. The fluid state is then extended out in the 
pitch direction by the equations given in the next section (compare with 
references 14 and 15). 


Variation of Fluid State in. Pitch Direction 


Equations (8), (10), (ll), and (12) directly give the first-order 
partial derivatives of W z , W , and p with respect to y in terms 

of the known quantities on the chosen (mean) streamline as follows: 



2 „ 
COS P 


cos 2 ^ 


+ w 


dW v 

y3T. 


(13) 

(14) 

(15) 


The second-order partial derivatives of W , W , and p with 

z j 

respect to y can he obtained as follows: Differentiating the con¬ 
tinuity equation (l) with respect to y results in 


d 2 (pw z ) 

dz dy 



(16) 


Equation (16) can be written through the use of relation (9) as 


d (pWy) 

dy 2 


o 2 (p«J d a (p'"V) 
- tan » ~~d/~ * to 


= 0 


(17) 


which is expanded to obtain 


\ , , 8p Sw 

2 Sy3T‘ 





tan p + — 
dz 


a( P w z ) 

~ 0 


» 


(18) 
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From equation (10) 


dW y _ . . 5w z 1 d (P W z) 

W~ ~ p 5y~ " p dT" 


Differentiating the irrotationality equation (2) with respect to y and 
using equation (9) yield 

*X d av y d 2 '" - ,, 

O- - 3Z ST - tan 13 v/ ( 2 °> 


Substituting equations (19) and (20) into equation (18) results in 


b \ 2 d (pW z ) Sp (d SW y _ b\\ d 

p dy 2 “ p Sr-Zy - p [teW ■ tan P g^~J tan P + di 

Transposing and combining terms give 


d(pWj 


2 d(pW Z } Sp 1 d f SW z ^ d aW yl 2 

2 ~ p as i p sr + w z ^ + tan P d^^r cos p 


2 bp 1 d 


After this equation is evaluated., the second partial derivative of W, 
with respect to y is obtained by using equation (20) : 


1 * ow 

z _ h y 


2“ " IS ' tan p 


A typical computation for these derivatives is presented in tables I 
and II. 


The second-order partial derivative of p with respect to y is 
obtained again from equation (12): 


C ^ B * W (h)(v? 


1 S 2 p_2-r/dpf l „ 


Sr; 


Third and higher order y-derivatives, if required, can be obtained 
in a similar manner. The complete variation of any fluid property q 
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across the channel can then be expressed by a Taylor's series in (y-y m ) 
from the various derivatives at the given streamline,, such as the mean 
streamline 



(24) 


Determination of Blade Profile 

The blade profile can be obtained by a consideration of mass flow. 

At the chosen stations, mass flow across a constant z line from y 

^m 

to y is computed as a function ,of y according to the following equa¬ 
tion: 


''y 

M = pW z dy (25) 

^m 

The variations of mass flow M and the magnitude of resultant velocity 
W at each station are plotted against y (fig. 2). Because the con¬ 
dition on the suction surface is more critical than that on the pressure 
surface, the blade shape on the suction side is determined first. From 

the plot of mass flow against y, a number of y T s are chosen for a 

s 

number of mass flows in the neighborhood of one-half the inlet mass 
flow, thus obtaining a number of suction surfaces. The corresponding 
velocities on the suction surfaces are read from the velocity plots 
(fig. 2). The one with the best velocity distribution is then chosen. 
After the suction surface is selected, the pressure surface and its 
velocity are determined in a similar manner by the total mass-flow 
requirement. If the shape or thickness of the blade or the velocity 
distribution on the blade obtained is not quite the one desired, the 
shape of and the flow on the mean streamline can be modified accordingly 
and the process repeated. Because each case takes a relatively small 
amount of computation, this modification is practical. Systematic 
building up of families of blades for various applications is also not 
difficult. 

The accuracy of the blade coordinates obtained depends mainly on 
the accuracy of series representation and the accuracy to which the 
partial derivations are evaluated. For high-solidity blades, such as 
those investigated in reference 13, the first three terms in the series 
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will give sufficient accuracy. For low-solidity blades, more terms may 
be required. It may be noted that accurate representation of the flow 
variation by the series is difficult to achieve in the neighborhood very 
close to the leading and trailing edges; but from a practical point of 
view it is satisfactory to fair in the nose and the tail according to 
some standard shape after the blade coordinates are obtained up to a 
short distance away from these regions. 

The application of this design method will be greatly aided by the 
availability of detailed flow variations in typical bladings such as 
those given in references 13 and 16. If a typical solution of the type 
of blading to be designed is not available, either a direct problem may 
be solved first, or even better, the accuracy of the inverse solution 
can be ascertained in the manner described in the following section. 
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Method of Checking Solution 


The inverse solution obtained by this method can be very conven¬ 
iently checked and improved, if necessary, by the relaxation method 
utilizing the fluid state obtained in the solution. Inasmuch as the 
velocities are available in the solution, the equation for irrotational 
absolute flow (equation (A9) of reference 13) is now written as 


dz 2 &y 2 


(«; 


8p 

y 5z 


W, 


dp 


= 0 


The finite-difference form of equation (26) is then 


(26) 



(v 1 o 3 

J=0 J 


(Wz ) 1 C 


fc =0 


1 i 
B 

n k 


k 

P 


0 


(27) 


where the same notation used in reference 13 is employed. A grid 
system is obtained by retaining the same z-stations used in the inverse 
solution and dividing the pitch distance into an appropriate number of 
divisions. The values of \J/, W, and p are most conveniently obtained 
by reading off the plots of integrated mass flow M, ¥, and p at 
each of the z-stations. The differentiation coefficients B's for 
equally spaced grid points as given in reference 17 can be applied 
throughout the domain for the present purpose by using the function 
values which are inside the blade but at equal spacing from the points 
in the channel (the first and last rows in tables III and IV). If the 
residuals obtained according to equation (27) require negligible change 
in \|/, the solution is entirely satisfactory. If the residuals are 
large enough to necessitate one cycle of Relaxation, the net effect 
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may be a slight change in the specified mean streamline flow and in the 
velocity distribution on the resultant blade. If the residuals are so 
large as to warrant a few cycles of relaxation, the flow variation for 
this type of blading is established, which makes the design of other 
similar bladings very simple. In general, for problems in which some 
knowledge of the flow is available, no relaxation should be necessary 
except, perhaps, near the nose and tail when accurate detailed velocity 
distribution in these regions is desired. 


Special Case of Incompressible Flow 

When the density is constant, the channel width ratio P/(P-t) has 
a relation to W z ^ m similar to the relation it has to (pW z ) m in the 

compressible case (reference 13). A number of values of W z are there¬ 
fore prescribed at a number of chosen stations along the mean streamline 
to lead to a desirable thickness distribution of the blade. The solu¬ 
tion of the incompressible problem continues in generally the same 
manner as it did for the compressible case with considerable simplifi¬ 
cation in the series terms and the integration process. In the incom¬ 
pressible case, of course, equations (12), (15), and (23) relating p 
and the velocities are unnecessary. 


The first- and second-order derivatives expressed by equations (13), 
(14), (21), and (22) for the compressible case are simplified to: 


and 


Sw y 

W 

dw z 

By" 


= \tan p 


aw, 

dz 


dW 

dz 


dW 

dz 

dW 


+ tan p 


dz 


z \ 2 

cos p 


z i 2 n 
COS p 


6% 


y 


By 


tan R ±(^l) . ±(^ 
^ dz \dy / dz \ 6y ‘ 


cos 2 p 


0 2 W 


z _ _d_ 
dz 


(SO- 


tan p 


Vw ; 

'By 2 


(28) 

(29) 

(30) 


(31) 


for the incompressible case. The equation for obtaining mass flow 
becomes M = p W z dy. 
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Illustrative Examples 

The procedure outlined in the preceding sections has been applied 
to the design of several highly cambered thick turbine blades for either 
compressible or incompressible flow. 

First and second examples. - In the first example, the shape of the 
mean streamline and the variation of axial velocity obtained in the 
incompressible solution of the blade given in reference 13 are taken 
as the specified values in order to determine whether the original 
-blade will be reproduced. The shape of the mean streamline is shown in 
figure 1 and the thickness distribution of the blade is shown in fig¬ 
ure 3. As a further check of the method, the z-stations chosen in this 
calculation correspond to every other z-station used in reference 13, 
so that the velocities obtained along these stations can also be com¬ 
pared with the solutions obtained in reference 13. The given quantities 
were inlet angle 0 ± , 41° 18'; exit angle 0 e , -52° 57'; axial chord 

L, 1.5 inches; and pitch P, 1.017 inches. 

The same data are used in the second example in which a compres¬ 
sible solution for an inlet Mach number of 0.42 is obtained. The axial- 
velocity variation on the mean streamline used in the first example is 
now taken as the specific-mass-flow variation on the mean streamline. 
This example is presented mainly to illustrate the difference in incom¬ 
pressible and compressible solutions for the same mean streamline shape 
and the same ratio of specific mass flow. 

The complete computation for the flow on the mean streamline and 
the determination of the first- and second-order derivatives of W , 

2j 

W , and p in the pitch direction at the mean streamline for these 

i/ 

two cases is given in tables I and II. Only three terms in the series 
are used because the direct solutions given in reference 13 indicate 
that they will be sufficient for engineering accuracy. The central- 
point fourth-degree differentiation formula is used at the regular 
stations z^ to z^. Because the first and last stations inside the 

channel employed in reference 13 are not close enough to the leading 
and trailing edges, respectively, two extra points are computed at 
stations 6.75 and 12.25 by the use of the unequal interval differen¬ 
tiation formula given in reference 18. These two points are so desig¬ 
nated because they are located at a distance of a quarter of the regular 
spacing away from stations 7 and 12, respectively. This combination is 
unnecessary for other cases. 

The velocities and the densities used in the calculation are non¬ 
dimens ionalized as follows: 
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W * = 


W 


w 


z,i 


W * = 


w_ 


w 


z,i 


* P 

° 


The mass flow M*, being divided by p.W ■, has the dimension of 
length: ' 


M* = 


M 


p.W . 
Z; 1 


2 

For the compressible case, the ratio of to W z ^ is equal 

to 25.78, which is the value- used in the construction of the two density 
graphs involved in the calculation, the reduced versions of which are 
shown in figures 4 and 5. 


The values in columns 2, 5, 8, 9, 12, and 13 of table I and 7, 8, 

13, 14, 22, and 23 of table II are used to compute the variation of W z 

and W in the pitch direction by the power series. The calculation 
«y 

of density follows directly as does the integration across the channel 
for mass flow, the mass flow being determined numerically. Because the 
specific mass flow pW_ was made dimensionless with its inlet value, 

the numerical value for the mass flow is equal to the pitch, the height 
being considered unity. This indicates that the integration for mass 
flow along the pitch direction starting at the mean streamline was con¬ 
sidered to be complete in either the plus or minus direction when a 
value equal to one-half the pitch distance or one-half the total mass 
flow was reached. Thus the channel flow and the blade coordinates 
were obtained. This calculation at one station, 10, is shown in 
tables III and TV for the two cases, respectively. 

The blade profile obtained by interpolating y g and for 

one-half of the inlet mass-flow value is shown in figure 6. Because of 
the same mean streamline shape and the same variation of specific mass 
flow on the mean streamline prescribed in the incompressible and com¬ 
pressible cases, the blades obtained for these two cases look quite 
similar except that the compressible one is somewhat thicker and the 
suction surface of the compressible solution is situated farther from 
the given mean streamline; these results are consistent with the direct 
solutions given in reference 13. The velocities obtained in the two 
cases (fig. 7), however, are quite different. Those in the compressible 
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solution are, in general, higher than those in the incompressible case, 
principally because of the high velocity on the mean streamline result¬ 
ing from the use of the same (pW.)^ and a decreasing p along the 

2j III 

mean streamline in the compressible case. 

The blade coordinates obtained in the incompressible solution are 
compared with those of the original blade in figure 8. The velocities 
at the six regular stations in the channel are compared (figs. 9 to 11) 
with the values obtained in the relaxation solution of the original 
blade reported in reference 13. These four figures show that, in 
general, the present solution is satisfactory. The relatively large 
differences near the leading and trailing edges are partly due to the 
inaccuracy in the second-degree polynomial approximation in the present 
calculation and partly due to the inaccuracy of the numerical solution 
obtained in reference 13 caused by the relatively coarse grid used in 
these regions. 

As a check of the consistency of the solution, an integration of 
the velocity along the blade profile obtained in the same example was 
made and compared with the circulation value computed from the inlet 
and exit tangential velocities and the pitch. The two are in agreement 
within 1 percent. 

Third and fourth examples. - In the previous two examples, the mean 
streamlines and specific-mass-flow distributions prescribed were not 
entirely arbitrary, having been obtained from the results of a direct 
problem for incompressible flow. The possibility of obtaining an 
unrealistic blade shape was, for this reason, largely eliminated. Con¬ 
sequently, in order to give the method a still more rigorous test, the 
thickness distribution (fig. 12) and mean blade line j(fig. 13) for the 
hub section of an experimental cooled turbine blade were arbitrarily 
chosen from which a mean streamline and specific mass flow were estima¬ 
ted by means discussed in the section "General Description of Method." 

In addition to the mean blade line and thickness distribution, the 
following data were used: inlet Mach number, 0.42; inlet angle, 36.2 ; 
exit angle, 42.7°; axial chord, 1.8 inches; and pitch, 1.176 inches. 

The mean blade line was faired into the inlet and exit directions 
(fig. 14) with modification to obtain a mean streamline according to 
the information obtained for a typical turbine blade in reference 13 
(see fig. 8). The thickness distribution was used to obtain the 
specific-mass-flow distribution along the mean streamline (dashed curve 
in fig. 15) according to results obtained in reference 13. After the 
flow on the mean streamline was calculated, it was then extended out 
from the mean streamline across the channel by means of power series as 
before. 

Both compressible and incompressible solutions were, completed, the 
incompressible case being considered first because of its relative 
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simplicity. The results of the incompressible case served as a guide to 
obtaining better prescribed values for the compressible case. The blade l 

obtained in the incompressible solution has a slightly higher thic kne ss 
distribution than was wanted and, because results obtained in the first 
two examples indicate that the compressible blade will be thicker than 
the incompressible for the same specific-mass-flow distribution, this 
distribution was depressed by a linear proportioning (fig. 15) in order 
to obtain better starting values for the compressible case and conse¬ 
quently a thinner blade (fig. 14). The same mean streamline numerically 
smoothed to give small fourth differences was used in both cases. 

As a check of the accuracy of the last solution, a grid having the 
same spacing in the z-direction used in the inverse solution (0.18) and 
a grid spacing of 0.147 in the y-direction is chosen, and the residuals 
at the grid points are computed according to equation (27) using the 
central point second-degree differentiation formula. As shown in 
table V, the residuals are rather small when they are compared with the 
magnitude of the coefficient at the points (-154.3). As an i nd ication 
of the percentage error in the velocity, these residuals are first 
divided by -154.3, resulting in an approximate change in the \|r value 
at each point. Then the probable error in dtyfd y or pW z is calcula¬ 
ted. The result is shown in table VI, which indicates that the solution 
obtained is sufficiently accurate for ordinary purposes. 


DESIGNING BLADES FOR COMPRESSIBLE FLOW ALONG ARBITRARY 
STREAM FILAMENT OF REVOLUTION 
Basic Relations 

The blade design method presented in the section "Designing Blades 
for Compressible Flow in Plane or on Cylindrical Surface" can be very 
easily extended to the more general case of flow along an arbitrary 
stream filament of revolution having a varying normal thickness 
(fig. 16). When only an average value in the stream filament of revolu¬ 
tion, as represented by the flow on the mean stream surface of revolu¬ 
tion described by a set of orthogonal coordinates 2 and cp (fig. 16), 
is considered, the equation of continuity for steady relative flow and 
the equation of irrotational absolute flow are given in reference 13 as 
follows: 


S(TpW z r) 

52 — 


5 (t pW^) 

59 


= 0 


dW„ A 

,je 
5 ~ 



W A sin o 

js. _ 


(i*) 
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+ 2o> sin a = 0 


( 2 >) 
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where primes in equation numbers indicate equations similar to those of 
the cylindrical case. 


When the fluid flow along a streamline S(l,cp) = 0 on the mean 
surface of revolution is considered, it is convenient to write any 
quantity q on S as a function of l only. Then the total deriva¬ 
tive of q with respect to l, following the streamline, is 


dq _ Sq L Sq dcp _ Sq , Sq 1 W cp _ Sq , tan 3 Sq 


dl ST Sep dl Si Sep r W 7 ST 


Sep 


(9') 


when 


¥ 


tan p 


= _ffi. = r 


w. 


dl 


When the preceding relation is used, the continuity and irrotationality 
conditions can be written as 


d(TpW z r) 

dl 


tan p 


dfrpW z ) dQrpWp) 


"Sep 


Sep 


= o 


( 10 ’) 


and 


dW 


cp 


d l 


r dcp 


tan 3 
r dcp 


W 


cp 


+ 2o> J sin a ~ 0 


( 11 ’) 


The variation of density throughout the flow region, in general, 
is obtained from the velocity by using the following equation: 



1 


A 12 2 

f 1 + 2 (av 

\ hi 



( 12 ') 


In order to obtain the density on the starting (mean) streamline, equa¬ 
tion (12’) is written in the following form (compare with equation (7a) 
of reference 13); 


2 



S> 


where 


(12a-) 
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and 


z 


*2 

P 



1 2 2 \ 
+ 2 mr } 
hi / 


2 

r-i 


® = ( P *w z * 


sec 




2 2' 
0) r 


h i 


r+i 
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Once the general relation between Z and $ is available, the evalua¬ 
tion of p* for certain given values of p*W^* sec (3 along the chosen 


(mean) 


streamline is made simple if auxiliary tables or graphs giving 

2 r+1 


/ T 1 2 2 \ ^-1 

(^) 

are first obtained, 


and 



I + 


2 2 
cd r 


h i 


r-i 


as functions of 


r 


The determination of the flow along the starting (mean) streamline 
(fig. 17) proceeds very much the same as in the case of plane flow or 
flow on a cylindrical surface. The shape of the streamline gives 
sec (3. The variation of p*W^* is obtained from the blade-thickness 
variation along the mean surface of revolution as follows: If the 

value of pW on the mean streamline represents its average value in 

L 

the circumferential direction 

T(pWj) m (P-t) = (TpW l P) i 

But 


P r 

hence 

T (p w z)m r = P 

Although there is always some deviation from this simple relation, 
especially around the leading and trailing edges, a relation between the 
two terms in equation (34) similar to that between (p*W *) m and 



(33) 


(34) 


(32) 
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P/(P-t) in the previous cylindrical case can be expected. Then, from 

T*(p*W *) r *. (p*W *) is calculated and combined with sec p from 
l m l m 

which is obtained by equation (12a 1 ). After p m * is determined, 

W, * and then * are easily calculated. 

Z,m cp,m 

Variation of Fluid State in Pitch Direction 

The first-order partial derivatives of and p with 

respect to cp are readily obtained from equations (10’), (ll*) > and 
(12*) as: 
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In a manner similar to that of the case of plane or cylindrical 
flow-the second-order partial derivatives of and p are 


obtained as follows: 
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After the variation of fluid state in the circumferential direction 
is determined, the mass flow across a constant l line from cp to 9 
is computed hy m 


M = rT 


tP 


<P, 


p W z dCp 


(25 1 ) 


m 


from -which the blade coordinates are determined as in the previous case. 
If checking is desired, the following equations are to be used: 
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where the same notation used in reference 13 is employed. 


SUMMARY OF RESULTS 

A rapid method for designing turbomachine blades of a given turning 
and a desirable blade-thickness distribution for a compressible non- 
viscous fluid flow along an arbitrary stream filament of revolution is 
presented. The method utilizes the guiding effects of the blade shape 
on the mean streamline shape and of the blade thickness on the specific 
mass flow along the mean streamline. After the flow on the mean stream¬ 
line is determined, the extension of the solution from the mean stream¬ 
line to the blade surfaces is accomplished by the use of a power series. 
A number of blade profiles are obtainable for the total mass-flow 
requirement, and one is chosen for the best velocity distribution on the 
blade. The results obtained in the solution can be used for a direct 
check on the accuracy of series approximation and, also, for the more 
accurate determination of the velocity distribution along the leading 
and trailing edges of the blade. 
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The method is illustrated with the design of several cascades of 
airfoils. In each case either some indication of the accuracy of the 
solution is given or the relation between the incompressible and com¬ 
pressible problems are shown. In the first case, the results obtained 
in an incompressible direct solution were used as prescribed values for 
the inverse problem to show that the original blade can be relatively 
easily reproduced with adequate accuracy. The variation of the velocity 
and its components across the channel was shown to compare favorably 
with those obtained in the direct solution and an integration around 
the blade for the circulation checked with that obtained from inlet and 
exit values within 1 percent. In the second problem, the same mean 
streamline and specific-mass-flow distribution along it were used; but 
this time the density was allowed to vary as in compressible flow and a 
comparison was made between the blades obtained in the first two cases. 

In the third and last problem, an arbitrary mean streamline and specific- 
mass-flow distribution were chosen, the blade obtained, and the solution 
checked by calculation of the residuals as in a direct relaxation solu¬ 
tion. Because the residuals were small the solution was considered to 
be reasonably accurate. A computation designed to give some indication 
of the accuracy of the velocities by dividing the derivative of the 
error in the stream function \|/ by the derivative of \|/ itself showed 
them to be accurate generally within 1 percent except near the blade 
boundaries where an end-point formula was necessarily used to obtain 
derivatives. 

Because the surface of revolution, on which the blades are located, 
is completely arbitrary, the method can be applied to axial-flow, radial- 
flow, and mixed-flow turbomachines. The variation in the normal 
distance between the stream surfaces of revolution can be taken into 
account, thus incorporating into the design the principal effect of 
three-dimensional flow. The method is readily applied to the design of 
channels on a plane and on a general surface of revolution. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, May 31, 1951. 
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TABLE I - COMPUTATION OF FLOW ON MEAN STREAMLINE AND FIRST- AND SECOND-ORDER y-DERIYATIVES AT THE MEAN STREAMLINE OF FIRST 

(INCOMPRESSIBLE FLOW) 
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TABLE II - COMPUTATION OF FLOW ON MEAN STREAMLINE AND FIRST- AND SECOND-ORDER 
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DERIVATIVES AT THE MEAN STREAMLINE OF SECOND EXAMPLE (COMPRESSIBLE 
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TABLE V - RESIDUALS OBTAINED IN COMPRESSIBLE SOLUTION OF THIRD EXAMPLE 



0.18 

0.36 

0.54 

0.72 

0.90 

1.08 

1.26 

1.44 

1.62 

1.80 

1.176 

1.029 

.882 

.735 

.588 

.441 

.294 

.147 

0 

- 1.2212 

- 1.0128 

- .3648 
.1086 

- .3708 

- .5373 

0.2093 
- .1610 

- .2084 
.2068 
.2642 

- .3869 

0.7716 
.3052 
.1641 
- .0093 
.0211 
.6299 

- 0.0361 

.3007 

.0091 

- .1896 

- .0233 

- .1013 

- 0.2726 

.1913 

.1671 

- .1197 
.2706 

- .1494 

0.0989 
.1127 
.0279 
- .1168 
.0586 
.0229 

0.0752 

- .3206 

- .2642 
.5123 

- .0189 

- .0003 

- 0.4316 
.2230 
- .1612 
.1244 
- .9797 
.2631 

- 0.0797 

- .2538 

- .1942 
.1099 
.0698 

- .2363 

- .4918 

- 0.2408 
.3271 
.4335 
.2336 
- .3507 
.0368 
.3491 


TABLE VI - ESTIMATED PERCENTAGE ERROR IN di|r/3y or pW z 



0.18 

0.36 

0.54 

0.72 

0.90 

1.08 

1.26 

1.44 

1.62 

1.80 

1.176 

1.029 

.882 

.735 

.588 

.441 

.294 

.147 

0 

- 0.00022 
.01599 
.01997 
- .00010 

- .00107 

- .00015 

- 0.01005 

- .00771 
.00644 
.00785 

- .00952 

- .01544 

- 0.01217 

- .01115 

- .00541 

- .00228 
.00941 
.01228 

0.01223 

.00081 

- .00830 

- .00051 
.00128 

- .00268 

0.01342 

.00779 

- .00524 
.00162 

- .00043 

- .01121 

0.00119 
- .00128 

- .00392 
.00049 
.00210 

- .00199 

0.01148 

- .00599 
.01381 
.00393 

- .00795 
.00437 

0.02231 

.00498 

- .00177 

- .01441 
.00236 
.04036 

- 0.00597 
- .00226 
.00691 
.00496 

- .00650 

- .01052 

- .00432 

0.01664 

.01366 

- .00190 

- .01598 

- .00394 
.01379 
.00555 
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Figure 2. - Determination of blade coordinate and velocity on blade. 
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Figure 5. - Variation of p* with W/W„ * in accordance with 

equation (12). 
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Figure 6. - Inverse solutions for compressible and incompressible flows. 
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Figure 8. - Comparison of incompressible solution -with 

original blade. 
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Figure 9. - Comparison of axial velocities. Vertical lines are stations. 
Using scale, read horizontal distance from station to corresponding 
curve to obtain velocity. 
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Figure 10. - Comparison of tangential velocities. Vertical lines are stations. 
Using scale, read horizontal distance from station to corresponding curve to 
obtain velocity. 
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Figure 11. - Comparison of resultant velocities. Vertical lines are stations. 
Using scale, read horizontal distance from station to corresponding curve to 
obtain velocity. 
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Figure 12. - Blade-thickness distribution. 
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Figure 13. - Desirable mean blade line and corresponding 

blade shape. 
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Figure 15. - Channel width ratio and specific mass flow along mean stream line used in incompressible and 

compressible solutions. 



































Figure 16. - Flow on arbitrary surface of revolution and stream filament of revolution 
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Figure 17. - Blade section on arbitrary surface of revolution. 
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